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The interaction of 32.5 and 6 nm ultrashort X-ray pulses with the solid materials B4C, SiC and 
Si is simulated with a non-local thermodynamic equilibrium (NLTE) radiation transfer code. We 
study the ionization dynamics as function of depth in the material, modifications of the opacity 
during irradiation and estimate crater depth. Furthermore, we compare the estimated crater depth 
with experimental data, for fiuences up to 2.2 J/cm 2 . Our results show that at 32.5 nm irradiation, 
the opacity changes with less than a factor of 2 for B4C and Si and a factor of 3 for SiC, for fiuences 
up to 200 J/cm 2 . At a laser wavelength of 6 nm, the model predicts a dramatic decrease in opacity 
due to the weak inverse bremsstrahlung, increasing the crater depth for high fiuences. 
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I. INTRODUCTION 

Novel laser sources that can generate extremely intense 
ultrashort X-ray pulses are developing fast. The Free 
Electron Laser in Hamburg (FLASH) [T] has produced 
soft X-ray pulses at lel4 W/cm 2 with a 20 micrometer 
focus at a wavelength of 32.5 nm, and will reach 6 nm 
by the end of 2007 [2J. At SLAC in California, the Linac 
Coherent Light Source (LCLS) [3] is expected to deliver 
ultrashort pulses of hard X-rays in 2008. 

Methods developed at large-scale facilities are likely to 
become more accessible due to the advance of table top 
sources with similar characteristics; high harmonic gener- 
ation [4, 5 has demonstrated 25 fs pulses at 10 14 W/cm 2 
[BJ, and laser wake-field accelerators producing mono- 
energetic electron beams [7] [5] IH] may act as a source 
for ultrashort X-rays. 

These lasers will access an unexplored regime for light- 
matter interaction, extending the capabilities of existing 
scientific techniques and enabling the development of new 
ones. Laser ablation is of great interest partly as a tool 
for processing material with nanometer precision, but 
also as a way of testing microscopic theories describing 
the interaction of X-ray laser pulses with solids. Predic- 
tive models could anticipate new applications and illumi- 
nate the limitations set by radiation damage in a sample 
and in focusing optics. Coherent diffraction imaging of 
biological specimens using free-electron lasers [TUJ, [TT] is 
an example of an application where the achievable res- 
olution depends on the ultrafast ionization dynamics of 
the X-ray-solid interaction. 

The interaction of ultrashort optical laser pulses with 
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solids has been studied extensively experimentally [12, 
[□] and theoretically [H E3 [16] . However, for VUV/X- 
ray pulses the dominant physical processes that govern 
the interaction are different, resulting in a quite different 
response of the material. Fajardo et al. [17 discuss the 
mechanisms of laser absorption for photon energies in the 
range of tens to a hundred eV, and pulse lengths of 1 ps 
to 1 ns. They present simulations of hydrodynamic ex- 
pansion for different pulse lengths and target materials. 
Theobald et al. [TH] has investigated the optical proper- 
ties of dense plasmas in the extreme ultraviolet spectral 
range using high harmonic generation. 

Recently, experiments at FLASH have been performed 
using a holographic time-delay technique |19j . The setup 
realized unique femtosecond time-resolved measurements 
of the structural changes induced by the FEL-material in- 
teraction. The modeling tools used to reproduce the ex- 
perimental results from [19l are presented in more depth 
by Hau-Riege et al. [2D], where the optical properties of 
dense plasmas are investigated assuming local thermody- 
namic equilibrium (LTE) conditions. 

In this article we take a step further and present non- 
local thermodynamic equilibrium (NLTE) simulations of 
the heating of solids by femtosecond VUV and X-ray 
lasers at different fiuences. Three different materials are 
considered: B4C, SiC and Si. These materials are in- 
teresting to compare due to their differences in Z and 
electronic structure, which result in quite different re- 
sponses to intense X-ray pulses. The materials are also 
candidates for X-ray optics, and the ablative properties 
of semiconductors are crucial in micro-processing appli- 
cations. First, we present simulations of VUV laser heat- 
ing experiments [2T] that were performed at FLASH. 
We define a threshold for ablation to identify ablation 
regions in the simulations and compare to experimen- 
tal crater depths. We then investigate how absorption 
mechanisms, and thereby depth profiles of heated tar- 
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gets, would change for the higher fluences which could be 
achieved at the FLASH-facility under tight focusing con- 
ditions. Finally, we present simulations for pulses around 
6 nm wavelength to investigate the ionization dynamics 
for higher photon energies. 



II. MODEL DESCRIPTION 

The interaction of intense ultrafast lasers with solid 
material is a complicated process. The plasma result- 
ing from such an interaction is a transient state of mat- 
ter that is exceptionally hard to probe. Ultrafast X- 
ray sources grant access to this regime by combining a 
short wavelength, which can propagate through the dense 
plasma, with a pulse length short enough to probe its 
temporal evolution [53] . Ultrafast X-ray sources are 
thus important both as primary pump pulses in applica- 
tions dealing with imaging [24] and material processing 
[2"5] . and as secondary pulses that can probe the region 
of interaction. A theoretical framework can either start 
from the solid state formulation (see e.g. [26 ) or from a 
plasma formulation originally developed for low-density 
hot matter [T7] [S7J . The former approach is suitable for 
fluences around the melting threshold, while the latter 
is suitable for higher fluence where ionization is consid- 
erable and chemical bonds have little significance. For 
the present study we choose a plasma formulation with 
high-density corrections. 

The laser-matter interaction is modeled in Cretin 
[28] 125] . a multi-dimensional NLTE radiation transfer 
code. The calculated level populations and transition 
rates provide opacities, heating rates and conduction co- 
efficients in each time step of the simulation. For this 
study we use the code's built-in screened hydrogenic 
atomic models, which are similar to those described in 
More 30 1 . The cold opacity has been scaled to fit the 
Henke absorption coefficients [31] at the laser wavelength 
for the specific materials. The electron energy distribu- 
tion is assumed to be thermal, i.e. the electrons thermal- 
ize instantaneously, so transition rates dependent on the 
electron energy distribution are obtained by integrating 
over a Maxwellian electron distribution. A study using 
the Fokker-Planck formalism [T3] for laser heated solids 
at similar temperatures and time scale concludes that the 
effect of a non-Maxwellian distribution is small. 

To account for the effects of high density, ionization 
potentials are lowered following the Stewart-Pyatt for- 
mula |32j . This is a commonly used approximate model 
and the results have been compared to those of more de- 
tailed models and to experimental emission spectra 33J . 
For solid densities, this continuum lowering has a consid- 
erable effect on the ionization dynamics. In the case of a 
32.5 nm, 25 fs laser pulse on SiC, the average ionization 
(zbar) is almost doubled at an intensity of lel5 W/cm2. 

The inverse bremsstrahlung process (free-free absorp- 
tion) applies the formalism of Dawson and Oberman |34j . 
The Coulomb logarithm uses a maximum impact param- 



eter b max = \J X D 2 + a>i 2 suggested by [35] for dense sys- 
tems. Xd is the Debye length and aj = (3/47rn i ) 1 / 3 is the 
ion sphere radius and n, is the ion density. The electron- 
ion coupling coefficient is calculated from Spitzer's for- 
mula [36,, using the same Coulomb logarithm as above. 

We model the heating in one dimension, justified by 
the laser spot size being much larger than the penetra- 
tion depth into the material. The simulations use 40 
zones which exchange energy through radiation transport 
and electron thermal conduction. The incident radiation 
is attenuated as it passes through each zone, and the 
opacity of each zone is calculated each time-step from 
the atomic populations and the inverse bremsstrahlung 
model. 

In the case of optical laser heating, a critical electron 
density for the propagation of the laser light is usually 
reached in an early phase of the the exposure and the 
laser energy is deposited mainly in a surface layer. The 
main channels of ionization are multi-photon ionization 
(MPI) and avalanche breakdown, where delocalized elec- 
trons gain energy through inverse bremsstrahlung fol- 
lowed by a geometric increase of free electrons through 
impact ionization. For intensities above m 10 15 W/cm 2 
the optical laser field becomes strong enough to suppress 
the atomic/ionic potential, causing field ionization. The 
Keldysh parameter (for a review, see e.g. [37]) indicates 
which process is dominant. For the short wavelengths 
considered in this study the Keldysh parameter 7 > 1, 
indicating that field ionization can be neglected. To esti- 
mate the degree of MPI, the simulations include below- 
threshold ionization with weak-field scalings from Delone 
and Krainov [3B"] . 

III. RESULTS AND DISCUSSION 

In this section we present comparison to experimental 
data for the materials B4C, SiC and Si, as well as pre- 
dictions of the response of the materials at higher fluence 
and photon energies than the current FEL beam param- 
eters. The theory presented in the modeling section is 
used to simulate the ionization dynamics throughout the 
material. The laser beam parameters have been chosen 
to match measured [T] and predicted [3pJ values at the 
FLASH facility. 

A. Crater depth and mechanisms of ablation 

Figure [I] shows the simulated average ionization (a) 
and electron temperature (b) as functions of depth for 
B4C at different fluences. The temporal structure of the 
pulse is chosen to match the pulse parameters at FLASH 
and is modeled by a Gaussian distribution with a FWHM 
of 25 fs, a wavelength of 32.5 nm, and with a spectral 
width of 1%. The intensity achieved so far in experi- 
ments at FLASH is similar to the red curves in figure [T] 
but we note that under tight focusing conditions, fluences 
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FIG. 1: The average ionization per atom (a) and electron 
temperature (b) after the laser pulse as a function sample 
depth for B4C at different fluences. 



in parity with the blue curve can be expected. The dam- 
age mechanism for VUV and X-ray irradiation in solid 
materials proceeds through different processes compared 
to that for optical lasers. In the latter case, valence elec- 
trons are excited to the conduction band, and as these 
free carriers accumulate they can be described as an elec- 
tron gas [26]. Lorentzo et al. [16] find that the optical 
laser-excited carriers at near-threshold ablation in silicon 
reach an equilibrium at around 0.8 eV through mutual 
collisions and impact ionization. On the other hand, for 
photon energies well above the binding energy of the va- 
lence electrons, photoelectrons are ejected into the lat- 
tice, generating many secondary electrons through im- 
pact ionization. This is a fast process, where an X-ray 
photon has generated a cascade of secondary electrons in 
a few femtoseconds [101 HI] • At a depth in the material 
where the beam is attenuated to a fluence approaching 
the threshold for ablation, the simulations suggest that 
the state of the material is similar to the optical case; the 
quasi- free electrons here, produced mainly through im- 
pact ionization, reach a temperature of 0.5-1.5 eV. The 
corresponding average ionization is around 0.1 electron 
per atom. The simulations also show that the electron- 
ion equilibration time is much longer than the pulse, ap- 
proximately 1 picosecond. Melting and ablation that oc- 



curs on a time-scale faster than conventional melting is 
usually referred to as non-thermal melting [26 . Calcula- 
tions [25] and measurements [32] for Si and GaAs indi- 
cate that the threshold can be expressed as the density 
of free carriers required to destabilize the lattice, and 
that its value is N u f « 10 22 cm~ 3 . In the simulations in 
this study, the depth where the density of free electrons 
reaches this threshold value corresponds to an electron 
temperature that lies within the range of critical tem- 
peratures of solids of 0.3-0.9 eV [43j [44] . This motivates 
us to identify this region of electron density and tempera- 
ture as a threshold for material removal. Thus, the simu- 
lations in figure [l] yield the depth into the material where 
we predict non-thermal ablation will occur. We note that 
recent studies of damage from the FLASH free-electron 
laser [25] [35] present experimental data supporting the 
idea that material removal proceeds via non-thermal ab- 
lation. 

Figure [2] shows a comparison between simulated and 
measured crater depths, with measured depths obtained 
from post mortem analysis of X-ray free-electron laser ir- 
radiated samples. The conditions of the experiment are 
described in more detail by Hau-Riege et al. [21] . The 
blue curves show that the predicted depths for B4C, SiC 
and Si. For B4C and SiC the agreement is good, even 
though the data is scarce above 1 J/cm 2 . For fluences 
close to the ablation threshold, the simulated depth is 
larger than the measured values for all materials. We 
note that the plasma formulation does not reproduce 
the thresholds for damage and ablation accurately; that 
would require a solid-state approach. In this study the 
aim is to quantify the response of the material during 
irradiation, and to estimate the crater depths for higher 
fluences. For Si the agreement is decent up to a fluence 
of 1 J/cm 2 , above which the experimental crater depths 
vary quite a lot. A similar effect was recently reported 
by Stojanovic et al. [35] , who observed clean craters in 
Si at 0.9 J/cm 2 , but a drastic increase in depth and a 
more complicated crater structure at 1.5 J/cm 2 . These 
findings, together with a time resolved study of optical 
reflectivity changes during heating and ablation, support 
the notion that the ejection of molten material becomes 
an important mechanism in Si above a certain fluence. 
They also point out that this "piston effect" is expected 
to be more pronounced for large absorption depths, which 
suggest that the such an effect would be less pronounced 
in B 4 C compared to Si. The crater depth data in figure 
|2jb) show a sudden increase and spread in depth above 
1 J/cm 2 . In line with the observations above, we assume 
that this effect is due to the expulsion of melted material. 
Hence, Si is excluded from the crater depth predictions 
at higher fluences. 

Figure [3] shows two quantities that are important for 
the ablation mechanisms in these materials. In (a) the 
deposited energy density at the end of a 1 J/cm 2 pulse 
is plotted as a function of depth. The steep profile for 
B4C (black curve) compared to Si (red curve) reflects the 
larger penetration depth of Si. The higher energy density 
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FIG. 2: Experimental crater depth as function of fluence (red) 
compared to heating simulations using the threshold for ul- 
trafast melting (UFM) (blue), for B 4 C (a), SiC (b) and Si (c). 
For B4C and SiC support the notion that material removal 
is initiated by ultrafast melting. In the case of Si, expulsion 
of melted material could be significant for fluences above 1 
J/cm 2 . 



at the surface in B4C is a result of the higher ionization 
and electron temperature. This leads to a hydrodynamic 
pressure that reaches 120 GPa, which is almost 6 times 
higher than that for Si (b). For both materials, the pres- 
sure at the end of the pulse is dominated by the electrons, 
with the ion and radiation pressure being less than 3% 
of the total pressure. The pressure at the B4C surface 
is higher, but for Si there is more material with a tem- 
perature around the melting/boiling point due to the flat 
heating profile. In combination with the lower heat ca- 
pacity of Si, this could possibly enhance the effect of melt 
expulsion. 
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FIG. 3: (a) The deposited energy density as function of depth 
for B4C (black curve) and Si (red curve) at the end of a 1 
J/cm 2 , 25 fs, 32.5 nm pulse. An exponential decay with depth 
indicates that the opacity of the material remains close to 
constant during the pulse, (b) Total electron pressure in the 
materials as function of depth. 



In the case of B4C and SiC, the amount of melt ex- 
pulsion appears to be considerably lower than for the Si 
case, and it has been suggested that the removal of ma- 
terial happens through two-phase vaporization followed 
by fluid-like expansion [21 . Large-scale hydrodynamic 
motion of 140 nm spherical polystyrene particles under 
the same beam conditions as in figure [2] has been con- 
firmed to take place during the first picoseconds |19j . 
During this time, our simulations indicate that further 
heating through electron thermal conduction and radia- 
tion transport from the hot region is negligible. However, 
such a sudden acceleration of material from the surface 
will release a compressive shock wave into the material 
[TTl |2T)1 |2"T] . To estimate the degree of heating in the ma- 
terial after the pulse, we simulate the material beneath 
the vaporized expanding material for 10 ps, and apply 
a Gaussian compression of the density (with a factor of 
2 as reported for similar conditions |17j ) over time as a 
crude model of a shock wave. The result suggests that 
there is very little further ionization after the pulse for 
the fluences considered in this study, and that the calcu- 
lated ionization depth-profile after the pulse can be used 
to estimate the crater depth when vaporization followed 
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by fluid-like expansion is the dominant mechanism. 



a) 



B. Predictions at high fluence 

From the analysis above, we see that the plasma model 
shows reasonable agreement with the experimental crater 
depths for B4C and SiC for fluences of 0.5-2 J/cm 2 . Si is 
included in the study of the material response during the 
pulse, but excluded in the prediction of the crater depth 
at high fluences due to the complicated ablative prop- 
erties. At higher fluences, the optical properties in the 
X-ray regime can be expected to change due to significant 
changes in the electronic structure. Such modifications 
can be studied by following the atomic populations in 
the constituent atoms. Below, we present predictions for 
the behavior of these materials during irradiation from 
pulses with fluences up to 500 J/cm 2 . This corresponds 
to the order of magnitude that can be expected at the 
FLASH facility under tight focusing conditions, and at 
future X-ray sources. 

The absorption of the laser light in the material is mod- 
ified due to plasma formation during the pulse. This is 
illustrated in figure [4] which shows the opacity at the 
laser wavelength of 32.5 nm as a function of depth at 
the end of an 200 J/cm 2 pulse for B 4 C (a) and SiC (b), 
and for a 500 J /cm 2 pulse for Si (c) . These fluences heat 
the materials to comparable electron temperatures at the 
end of the pulse (~ 125 eV). At a laser wavelength of 
32.5 nm there is no contribution from the bound-bound 
transitions for these elements, so the opacity comes from 
bound-free (bf) and free-free (ff) transitions. The right 
hand side of the plots corresponds to a depth where the 
material is cold and the opacity is proportional to the 
initial photoionization cross section. As we move left 
toward the surface, the bf-contribution decreases as the 
atoms ionize, and the ff-contribution increases as a result 
of the increase in the density of quasi-free electrons. In 
the case of B 4 C, only the core electrons remain in the 
zones closest to the surface, leading to a flattening of 
the opacity curve (bf+ff) with depth. The depth pro- 
file is similar for SiC. In the case of pure Si, the opacity 
starts to rise close to the surface. This is because the 
L-shell electrons in Si have a lower ionization potential 
than the K-shell electrons in C and B, producing a higher 
density of quasi-free electrons for the same temperature, 
and hence a stronger ff-contribution. It is of interest to 
find the fluence where the opacity starts to deviate from 
the cold value. For B4C and SiC the opacity changes 
by 25% at about 5 J/cm 2 . For Si the fluence is about 
20 J/cm 2 . We also note that the free- free contribution 
becomes larger than the bound-free contribution at elec- 
tron temperatures of 35, 31 and 23 eV for B4C, SiC and 
Si respectively. 

To emphasize the difference in ionization dynamics 
compared to an LTE treatment, figure [3] also shows the 
opacity resulting from a Saha-Boltzmann population dis- 
tribution at the same temperatures as in the non-LTE 
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FIG. 4: Opacity as function of sample depth at the end of a 
25 fs, 32.5 nm pulse. The fluence is 200 J/cm 2 in the case 
of B4C (a) and SiC (b), and 500 J/cm 2 in the case of Si 
(c). The non-LTE total opacity (bf+ff, red curve) is the sum 
of the free-free (ff) contribution (blue) and the bound-free 
contribution (black). The bf opacity curve converges to the 
cold opacity at large depths, whereas the ff opacity is sensitive 
to the number of free electrons and hence increase closer to 
the surface. For comparison, the dotted curve shows the total 
opacity assuming LTE, at the same temperatures as in the 
non-LTE case. The LTE opacity model results in a more 
restrained ionization at low temperatures, whereas at high 
temperatures the LTE ionization is stronger, increasing the 
inverse bremsstrahlung absorption. 



case. The main reason for the differences is that the 
LTE treatment results in a higher average ionization, 
leading to enhanced inverse bremsstrahlung absorption. 
The overestimation of the average ionization in LTE is 
known from previous studies |46l I47j . By performing the 
same analysis at lower fluences we can find where the 
departure from LTE begins to have an impact on the 
opacity at the surface. We choose to consider a differ- 
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ence in opacity of 25% as "possible to distinguish" in, 
e.g. a transmission measurement. This is reached at an 
electron temperature of 45 eV and 92 eV for Si and B 4 C 
respectively, validating the LTE treatment used in 20J. 
These temperatures are reached for similar fluences; 65 
and 90 J/cm 2 respectively. We note that this result holds 
for similar pulse lengths and wavelengths only; an esti- 
mate of the deviation from LTE must be found whenever 
the pulse parameters change significantly. 
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FIG. 5: The opacity (in cm 2 /g) for B4C (a) and Si (b) as a 
function of material depth and time during pulse irradiation, 
for the same pulse parameters as in figure [4] The materials 
reach a similar electron temperature at the surface at the end 
of the pulse. The opacity is rapidly reduced by about a factor 
of two in the surface region for both materials, but in the case 
of Si it reaches the initial value at the end of the pulse due to 
the higher rate of inverse bremsstrahlung heating. 

The heating depth profile after the pulse depends on 
the opacity throughout the interaction. Figure [5] illus- 
trates the temporal evolution of the opacity during the 
pulse as function of depth for B4C (a) and Si (b), with 
the same fluences as in figure [4] The SiC depth-time 
profile is similar to that of B4C. The material in the 
blue-colored regions attenuates less radiation, but due 
to the inverse bremsstrahlung contribution, the opacity 
doesn't decrease by more than a factor of two. In the 
case of Si, the density of free electrons is higher, even 
though the temperature in the region close to the surface 
is similar. This leads to increased inverse bremsstrahlung 
absorption, increasing the opacity at the surface during 
the last 10 fs of the simulation. Figure [5] also shows the 
depth at which a significant modification of the opacity 
can be expected. We note that by choosing a sample with 



a thickness similar to this depth (e.g. a thin film), one 
could measure the transmission of either the pump-beam 
or of a probe beam arriving at a later time. Such a set- 
up is straightforward in principle, and could be useful for 
testing NLTE predictions, or for validating existing mod- 
els for inverse bremsstrahlung in solid density plasmas. 
Furthermore, by probing the region of interaction after 
heating it with a second pulse, one could measure the 
electron density. Such a measurement could contribute 
to the process of validating hot dense plasma models. 
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FIG. 6: Estimated crater depth for B4C (red) and SiC (blue) 
based on the heating profile after the pulse, using the de- 
fined threshold for non-thermal ablation. The laser wave- 
length is 32.5 nm and the pulse length is 25 fs. The decrease 
in photoionization is almost compensated for by the inverse 
bremsstrahlung, resulting in an exponential dependence be- 
tween depth and fluence up to about 10 J/cm 2 . For higher 
fluences the crater depths increase faster due to the decreased 
opacity. 

Figure [6] shows the estimated crater depth for B4C 
(red) and SiC (blue) as functions of fluence based on the 
simulated heating profiles, using the threshold for non- 
thermal ablation. The region below 3 J/cm 2 is the same 
as the blue curves in figure[2]ja) and (b). The effect of the 
reduced opacity has a marginal effect on the estimated 
crater depth for fluences below 10 J/cm 2 , as confirmed 
by the exponential dependence on fluence. Assuming a 
constant opacity k, we expect an attenuation of the in- 
cident photon flux I : I = I e~ Kpd , where d the depth 
into the material and p is the density of the material. 
Such a dependence would result in a straight line for the 
log-plot of figure [6j We conclude that the rate of in- 
verse bremsstrahlung absorption is almost high enough 
to compensate for the decrease in the photoionization 
cross section. However, a moderate decrease in the total 
opacity can be identified above 10 J/cm 2 , especially for 
SiC (blue). 

C. Predictions at 6 nm wavelength 

In this section we compare the results of simulations in 
the vacuum-ultraviolet region with those of simulations 
in the mid-soft X-ray region. We choose a wavelength 
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close to 6 nm since this will be available for experiments 
in the next stage of the FLASH soft X-ray free-electron 
laser [T] . This study does not consider resonance heating 
where the laser wavelength overlaps with the absorption- 
peaks resulting from the bound-bound contributions. We 
hence choose to simulate the interaction of 5.9 nm radi- 
ation with B4C. 
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FIG. 7: Temporal evolution of the electron temperature at 
the surface with a Gaussian pulse profile (red) and a spiky 
pulse profile (black), at a wavelength of 5.9 nm, a flucncc 
of 200J/cm 2 and a pulse length of 15 fs. The correspond- 
ing pulse shapes are shown in blue (Gaussian) and magenta 
(spiky), with an arbitrary intensity. Both pulses contain the 
same number of photons, and result in the same final electron 
temperature. 

The temporal profile of the FEL pulse is expected to 
contain more pronounced spikes for shorter wavelength 
[59"] . To check how such a spiky profile affects the ioniza- 
tion dynamics, we simulate the heating of B4C using both 
a Gaussian pulse and a pulse with a spiky temporal pro- 
file made to resemble the predictions by Saldin et al. [39] ■ 
The two pulses have equal energy (200 J/cm 2 ) and pulse 
length (15 fs), and their shapes are shown in figure [7} The 
Gaussian pulse (blue) induces a time-dependent electron 
temperature described by the red curve, while the spiky 
pulse (magenta) corresponds to the black temperature- 
curve. The temperatures follow different paths during 
the first part of the simulation, but converge to the same 
temperature during the second part of the simulation. 
Hence figure [7] supports the use of Gaussian pulse as a 
valid approximation in the following analysis. 

Figure [8] shows the temporal evolution of the opacity 
as a function of material depth for a 5.9 nm, 200 J/cm 2 
Gaussian pulse. At this wavelength the photons are en- 
ergetic enough to ionize the boron K-shell directly, but 
not the carbon K-shell. The incident photons interact 
mainly with the core electrons in boron, with photoion- 
izations followed by fast impact ionization of the valence 
electrons of both boron and carbon. The free electrons 
are heated to a temperature of about 40 eV, which is too 
low for significant impact excitation or ionization of the 
K-shell electrons. In the case of 32.5 nm irradiation at 
the same fluence (described above), the photons interact 
with the valence electrons. The free electrons are heated 
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FIG. 8: Temporal evolution of the opacity as function of ma- 
terial depth for B4C. The pulse has a wavelength of 5.9 nm, a 
fluence of 200 J/cm 2 and a pulse length of 15 fs. The opacity 
decreases by almost a factor of four due to the weak inverse 
bremsstrahlung. 



to a temperature of about 100 eV. This temperature is 
high enough for the electrons in the high-energy tail of 
the electron energy distribution to cause some excitation 
and ionization of the K-shell electrons in both boron and 
carbon. The cold opacity at 5.9 nm is very similar to that 
at 32.5 nm, but during the pulse the opacity decreases 
by almost a factor of four, mainly due to the weaker in- 
verse bremsstrahlung absorption at this relatively short 
wavelength. 
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FIG. 9: Estimated crater depth for B4C based on the heating 
profile after the pulse, using the defined threshold for non- 
thermal ablation. The 32.5 nm case is shown for comparison. 
There is a drastic decrease in opacity due to depopulation of 
the boron K-shell and the weak inverse bremsstrahlung, lead- 
ing to an increase in the predicted rate of ablated material. 

The consequence of the decrease of the opacity at 5.9 
nm on the estimated crater depth is shown in figure [9] 
For fluences up to about 2 J/cm 2 the estimated heat- 
ing depths are very similar to the 32.5 nm case, but for 
higher fluences the rate of ablated material is predicted 
to increase. 

We note that in the resonant region of the absorption 
spectrum, there will be line contributions from bound- 
bound transitions in the plasma, further complicating the 
absorption mechanisms. We expect considerable broad- 
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ening of the lines at solid density, especially for low Z ma- 
terials like B4C. In the case of Si, we expect the lines to 
result in wavelength regions of significantly increased ab- 
sorption. Such an increase could affect the mechanisms 
of ablation and possibly increase the threshold flucncc 
for melt expulsion, allowing more control over the crater 
depth and providing a cleaner ablation at higher fluences 
(due to the suppression of melt expulsion) . 

IV. CONCLUSIONS 

The advent of intense femtosecond pulses in the X-ray 
regime opens up a new branch of laser physics, where the 
interaction with materials proceeds by different physical 
processes compared to the optical regime. Beside the in- 
herent interest of this regime, a quantitative understand- 
ing of the interaction is critical to find the optimal work- 
ing parameters for demonstrated potential applications 
[2"4"l |2"5] . We have addressed this issue by performing sim- 
ulations of the heating using a plasma code that follows 
the atomic populations in time throughout the material. 

In conclusion, we have simulated the interaction of 
intense ultrafast X-ray pulses with the materials boron 
carbide (B4C), silicon carbide (SiC) and silicon using a 
non-local thermodynamic equilibrium (NLTE) radiation 
transfer code. We chose pulse parameters to match recent 
experiments at FLASH [lj, and produced predictions for 



the next stage of operation of that facility [39 . Our simu- 
lations show that at 32.5 nm wavelength, the decrease in 
opacity due to photoionization is partly compensated for 
by increased inverse bremsstrahlung absorption, leading 
to a near-exponential energy deposition in the material. 
At a wavelength of 6 nm this effect weakens, resulting 
in deeper energy deposition for fluences above 2 J/cm 2 . 
Crater depths can be estimated from the simulated heat- 
ing profiles and the thresholds for ultrafast melting. To 
evaluate the model and investigate the mechanisms of 
ablation, we include a comparison to experimental crater 
depth data [ST] for fluences up to 2.2 J/cm 2 . The ex- 
perimental data for Si, and the ID heating simulations 
support the idea proposed by Stojanovic et al. [45] that 
melt expulsion leads to an increase in the crater depth at 
fluences above « 1 J/cm 2 . Our simulations also indicate 
that, for these pulse parameters, an NLTE treatment is 
required for fluences above approximately 80 J/cm 2 . 
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